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Abstract 

A theory is set up of spherical proteins interacting by screened electrostatics and constant adhe- 
sion, in which the effective adhesion parameter is optimized by a variational principle for the free 
energy. An analytical approach to the second virial coefficient is first outlined by balancing the 
repulsive electrostatics against part of the bare adhesion. A theory similar in spirit is developed 
at nonzero concentrations by assuming an appropriate Baxter model as the reference state. The 
first-order term in a functional expansion of the free energy is set equal to zero which determines 
the effective adhesion as a function of salt and protein concentrations. The resulting theory is 
shown to have fairly good predictive power for the ionic-strength dependence of both the second 
virial coefficient and the osmotic pressure or compressibility of lysozyme up to about 0.2 volume 
fraction. 
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I. INTRODUCTION 



ution properties of globular proteins may bear relation 
Since the characterization of proteins commands 



It has been intimated that the so 
with their crystallization properties 
ever more attention, such a contention is of considerable interest so much work has been 
carried out on this topic recently 

The difficulty of setting up a predictive theory of protein suspensions based on what is 
mown about the interaction between two proteins, has been acknowledged for some time 
jj. Best fitting of the osmotic pressure of, for instance, bovine serum albumin up to 100 
g / 1, leads to effective excluded volumes whose behavior as a function of salt is enigmatic 

In recent years, there has been a tendency to forget about all detail of the protein inter- 
action altogether — both attractive and repulsive — and simply introduce a single adhesion 
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Despite the electrostatic repulsion which is substantial, the 
data are often merely rationalized in terms of the bare protein diameter within 



parameter |l( 



die context 
12. 



of an adhesive sphere model and such an approach seems to have merit [10J, |11|, |12j, [13|, [l^ . 
This empiricism has prompted us to develop a theory of screened charged protein spheres 
that have a constant stickiness, where the electrostatic interaction is compensated, in part, 
by the adhesive forces. Thus, we argue that, effectively, the spheres are asigned a hard diam- 
eter identical to the actual diameter provided the remnant adhesive interaction now depends 
on the electrolyte and protein concentrations in a manner to be determined variationally. 
First, we analyze the second virial coefficient as such, for this will point toward a way of 
dealing with the osmotic pressure at nonzero concentrations. We focus on experiments with 
ysozyme, a protein which is reasonably spherical and has been well studied for a long time 



2 



II. SECOND VIRIAL COEFFICIENT 



A. Theory 

1. Second virial coefficient 

The second virial coefficient B 2 describes the first order correction to Van 't Hoff 's law 

-^- = l + B 2P + 0(p 2 ). (1) 
pk B T 

Here, II is the osmotic pressure of the solution, p is the particle number density, k B is 
Boltzmann's constant and T is the temperature. From statistical mechanics we know that, 
given the potential of mean force U (r) between two spherical particles whose centers of mass 
are separated by the vector r, one can calculate B 2 from 

B 2 = - l -j^dvf(v) (2) 

where /(r) = e ~ u ( r )/ k BT — 1 is the Mayer function. In principle, the interaction U(r) may 
be determined from experimental data on the second virial coefficient by suitable Laplace 
inversion. This has been done for atoms and spherically symmetric molecules Q, \v\ . 
for which the second virial coefficient has been measured over a broad enough range of 
temperatures. One might think of formulating a procedure similar in spirit and applicable 
to protein solutions, but with the ionic strength as independent variable instead of the 
temperature. However, to be able to determine an accurate approximation of the interaction, 
the experimental data have to be known fairly accurately, which is not the case at hand, 
as will become clear further on. We are therefore forced to adduce presumptions about the 
interaction. 

We assume the protein to be spherical with radius a, its charge being distributed uniformly 
on its surface. For convenience, all distances will be scaled by the radius a of the sphere and 
all energies will be in units of k B T. Because monovalent ions (counterions and salt ions) are 
also present in solution, there will be a screened Coulomb repulsion between the proteins, 
here given by a far- field Debye-Hiickel potential. We compute the effective charge qZ e ff in 
the Poisson-Boltzmann approximation where q is the elementary charge. For now, we let 
the attraction between two proteins be of short range, and we model it by a potential well 
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of depth U a and width S <C 1. The total interaction U(x) between two proteins is of the 
form 



U(x) 



oo < x < 2 

U DH {x) -U A 2<x<2 + 5, (3) 
Udh(x) x>2 + 5 



r 

x = — 







a 



U DH (x) = It- . (4) 



with Debye-Hiickel potential 

X 

Here, £ = ^ f-fq^H > K 1 is the Debye length defined by k 2 = 8nQI, I is the ionic strength, 
Q = q 2 /tkBT is the Bjerrum length, which equals 0.71 nm in water at 298 K, e is the 
permittivity of water and /z = kcl = 3.28a^/l in water, if a is given in nm and / in mol /\. 
We suppose 1-1 electrolyte has been added in excess so / is the salt concentration. 

In order to evaluate B 2 analytically, we have found it expedient to split up B 2 into several 
terms: 

B, = B^ s (l + |j) , (5) 

where B^ = 167ra 3 /3 is the second virial coefficient if the proteins were merely hard spheres 
and we introduce integrals 

POO 

J = J dxx 2 (l- e~ u{x) ) = J x - (e UA - 1) J 2 , (6) 

POO 

J x = / dx x 2 (1 - e~ UDH ^) , (7) 



2+<5 

J 2 = I dxx 2 e~ UDH{x) . (8) 



2 

Here, J\ is the value of J in the absence of attraction and may be simplified by Taylor 
expanding the Boltzmann factor in the integrand for small values of Uoh to second order. 
However, to increase the accuracy of the expansion, we adjust the coefficient of the second 
order term so that the approximation to the integrand coincides with its actual value at x = 
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2, i.e., we approximate x (l — e ^hM) ~ 2£e M( ^ x ' 2 ) — 2ct£ 2 e 2 ^ x 2 \ with a 
resulting in 



-(i-O 



(9) 



where we have neglected the small term a£ 2 /2/z 2 . For instance, in the case of lysozyme, the 
deviation of the approximation Eq. © from the exact result is smaller than about 3% for 
/ > 0.05 M and smaller than about 1% for / > 0.2 M. Since 8 <^1, J 2 may be simplified by 
using the trigonometrical approximation J^ +5 dx g(x) ~ \8 [g{2) + g(2 + 5)}, which leads to 



Jo ~ 25 



8 



+ 1 + - e T + 7 ^ 



(10) 



It is important to note that fi8 may be greater than unity even if 5 <C 1. Again, for lysozyme, 
this approximation deviates less than about 3% from the exact value for I > 0.2 M and 
5 < 0.5 and less than about 1% for I > 0.2 M and 8 < 0.15. 



2. Effective attractive well 

We next present a discussion of B 2 in terms of equivalent interactions and their Mayer 
functions even though the analysis of the previous section is self-contained. Sections IIA2 
and IIA3 may be viewed as preludes to the formulation of the liquid-state theory developed 
in section III. At large separations (x > 2 + 5), the interaction between the particles is purely 
repulsive, leading to a positive contribution to the second virial coefficient. If, at a certain 
ionic strength, the second virial coefficient is smaller than the hard-core value {B 2 < B^), 
this positive contribution is necessarily cancelled by only part of the negative contribution 
of the attractive interaction at small separations, the part, say, between x = 2 + €q and 
x = 2 + 8, see Fig. 1. The remaining potential which we will call an effective attractive well, 
then consists of a hard-core repulsion plus a short-range attraction of range eo- The value 
of 6q is determined by noting that the free energy of the suspension must remain invariant, 
which, in the asymptotic limit of low densities, leads to the identity 

B 2 ,e = B 2 , (11) 

where B 2 is the second virial coefficient of the previous section and B 2 ^ Q is the second virial 
coefficient pertaining to the effective attractive well. Using Eq. (J2J), we rewrite Eq. (fTTj) as 

/ d 3 r Af = 0, (12) 
Jv 
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in terms of the difference in the respective Mayer functions 



A/ = /-/ eo , 



(13) 



where / is the Mayer function of the original interaction and f eo is the Mayer function of 
the effective attractive well. In dimensionless units, Eq. (|T2j) is equivalent to the condition 



dx x 2 (l 



-U DH {x) 



'2+5 J2+e 

where, using the same approximation that led to Eq 



) = [ 2+S dx x 2 (e UA e- u ^ - 1) 

J2+en 



we write 



dx x 2 (l 



-U D h{x) 



) 



'2+5 ^ 

and, using dx x 2 Af(x) ~ 2(5 - e ) [A/(2 + 5) + A/(2 + e )] 



(14) 



(15) 



2+8 



2+e 



dx x 2 ( e UA e- UDH{x) - 1) ~ 2(5 - e ) 



2 + ( e 1+i / 2e ' + e 1+e o/ 2 



To leading order, we then find an explicit relation for So 

8-e 



jje 



u A 



(16) 



(17) 



which works well at high ionic strengths (i.e. at low values of £), e.g. whenever I > 1 M 
in the case of lysozyme at pH 4.5. A more accurate value of 5 — €q is obtained by equating 
Eqs. (|T5|) and (JTHJ), and then iteratively updating the factor (8 — e ), starting with the initial 
value e = 5. 

The second virial coefficient pertaining to the original potential U(x) (Eq. (J3J) is now 
rewritten as 

-2+eo 



B 2 = B 2 , = B? s (l + lj + " ^ x 2 (1 - e u -e- u °»^ . 



(18) 



The depth Ua — Udh(x) does not vary strongly though, since e$ <C 1, so, to simplify things 
computationally, let us approximate the interaction by a square well potential, 



oo < x < 2 
Usw(x) = { -U s 2 < x < 2 + e 
x > 2 + e 

We choose t/g in such a way that _E?2 = B% or > equivalently, 

"2+e 



(19) 



dx x 2 (e Us - e UA e -u DH (x)j = Q 



(20) 



To leading order in eo, we have 



J dxx 2 e Us ~4e e as , (21) 



and, using the approximation J 2 60 dx x 2 g(x) ~ 2e [g{2 + e ) + g(2)], we write 



x 2 e UA e~ UDH( ^ ~ 2e e^ 
The depth Us of the potential is then given by 



e 5 _)_ e i+ eo /2 e 



(22) 



~ L b Ua ( e -€ + e -i472 e - M£0 ) (23) 

in terms of the original variables. Finally, we point out that the two attractive wells that 
we have introduced are physically meaningful only if B 2 < B 2 S ■ 



3. Attractive well in the Baxter limit 

We have shown that one may simplify the statistical thermodynamics of the protein 
suspension at low densities considerably, by replacing the original interaction, consisting of 
an electrostatic repulsion and a short-range attraction, by a single attractive well of short 
range. The electrostatic interaction may be substantial but it is compensated by part of the 
original attractive well which is quite strong (Ua > 1)- Another useful interaction expressing 
attractive forces of short range consists of a hard-sphere repulsion and an attraction of infinite 



strength and zero range, namely the adhesive hard sphere (AHS) potential of Baxter 
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U A hs{x) = < 



oo < x < 2 

logg^ 2<a;<2 + ^ , (24) 
x > 2 + oj 

where r is a constant and the limit u J. has to be taken after formal integrations. The 
second virial coefficient remains finite 

B** S = B? s (l - -L) . (25) 

Because much is known about the statistical mechanics of the Baxter model, one often defines 
r in terms of some B 2 and naively assumes there is a one-to-one correspondence between the 
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original and Baxter models. For instance, in our case, b£ hs = b 2 = b 2 
we have 



Bl w . Since 



B» s 



B. 



HS 



1) 



1 



2/ 



- 1 



(26) 



l)e 



we thus identify 



6e (e Us - 1) , 



(27) 



where Us is given by Eq. ([23)1 . However, it is important to realize that this procedure is 
legitimate at small densities only. At finite concentrations, the optimal representation of the 
real suspension of proteins by a Baxter model has to be derived and we will show in section 
III that the simple-minded identification B^ 113 = B 2 no longer applies. 



B. Application to lysozyme 
1. Experimental Data 

Lysozyme is, by far, the best studied protein with regard to solution properties. This is 
one of the reasons for using this protein to test theory, another being its moderate aspect ratio 
of about 1.5 so that it may be fairly well approximated by a sphere. Bovine Serum Albumin 
(BSA) has also been well studied, but is considerably more anisometric with an aspect ratio 
of about 3.5. Numerous measurements of the second virial coefficient of lysozyme have 



)een publis 



14, 



20 
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red. In fact, there are quite a few sets of experiments pertinent to our analysis 
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24, 



25 



26 



21M 



It turns out that there is appreciable scatter in the data if we plot all measurements of 
B 2 at a pH of about 4.5 as a function of ionic strength / (NaCl + small amount of Na 
acetate; we have set the ionic strength arising from the latter equal to 0.6 x concentration 



2l|) (see Fig. 2). Several sets of data 



25 



281 ] appear to be way off the general curve within 



any reasonable margin of error. An important criterion is how well the 9 point (i.e. when 
B 2 = 0) is established since then attractive forces — which we would like to understand — are 
well balanced against electrostatics — which we purportedly understand well. Experimentally 
speaking, it ought to be possible to monitor B 2 accurately about the 9 point; large negative 
B 2 values at I ^> I g are more difficult to determine because the proteins may start to 
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aggregate or nucleate, in principle. Various polynomial fits for a 



1 data close to the 9 point 



yield Iq = 0.20 ± 0.01 M. Hence, we have regarded data sets [2a, |28[ markedly disagreeing 
with this ionic strength as anomalous so we have not taken them into consideration. Fig. 
3 displays all data we have taken into account. Clearly, the composite curve yields a fairly 
reliable basis to test possible theories of the attractive force. On the other hand, it is unclear 
at present how the scatter in data in Fig. 2 translates into bounds for inferred attractive 
interactions. 



2. Theory 

2a. Electrostatics 



Next, it is important to ascertain the actual and effective charges of lysozyme under 
conditions relevant to the present work. Kuehner et al 2^ performed hydrogen-ion titra- 
tions on hen-egg-white lysozyme in KC1 solutions. By interpolation, we obtain the actual 
charge Z of the protein as a function of the 1-1 electrolyte concentration / (see tables I and 
II). Experiments on B 2 are usually carried out with NaCl (and some Na acetate) as the 
supporting monovalent electrolyte but we here assume KC1 and NaCl behave identically in 
an electrostatic sense. We solve the Poisson-Boltzmann equation to get the effective charge 
Z e ff in the Debye-Hiickel tail (for more detail, see Appendix A). The dimensionless radius 
is set equal to \i = 3.28a-\/7 = 5.58-\/7 and Eq. ()61|) is used to compute the renormalized or 
effective charge. (Setting a = 1.7 nm for lysozyme as in Refs. (20] and j^j]; the Bjerrum 
length Q = 0.71 nm for H2O at room temperature). The other dimensionless parameter is 
given by f = 0.209(^/(1 + fj,)) 2 , where Z = Z eff - 1 (see below). 



2b. Attractive well 



We have assumed Ua and 5 to be independent of the ionic strength I. It is possible to 
show that this does not contradict the data displayed in Figs. 2 and 4. In Appendix B, we 
prove that if the interaction between the proteins is given by Eq. (JHJ), then dB2/dfi < and 
d 2 B 2 /dfi 2 > 0, the last inequality being valid if £ < 1. We recall that fi is proportional to 
\fl so that Figs. 3 and 4 indeed bear out these inequalities after due rearrangement. 
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Next, we determine the optimal values of Ua and 5 yielding exact, numerical -82(7) 
curves given by Eq. (jSJ) which are the best fits to the data of Fig. 3. We require that 
Iq = 0.20 ± 0.01 is predicted absolutely which fixes Ua, say, and 5 is then determined by a 
nonlinear minimization procedure. We thus obtain Ua = 1-70 ± 0.25 and 5 = 0.468 =F 0.097 
but we note that the quantity 5 exp Ua = 2.56 ±0.10 is much more narrowly bounded. Now, 
it can be argued that the Debye-Huckel potential with effective charge Z e ff overestimates 
the real potential in magnitude so we have repeated this numerical procedure with a slightly 
lower effective charge, viz. Z = Z e ff — 1 (see tables I and II). This yields the revised 
estimates U A = 2.87 ± 0.65, 5 = 0.167 =f 0.086 and 5 exp U A = 2.95 ± 0.21. The numerically 
computed curves are displayed in Fig. 3. We therefore conclude that the variables Ua and 
5 as such are difficult to ascertain unambiguously, though the variable 5 exp Ua is quite 
robust. This is also borne out if we use our approximations, Eqs. © and (JTUJ), instead of 
the exact numerical computations. There are again wide variations in Ua and 5 but the 
quantity 5expf/^ is strictly bounded: SexpUA = 2.70 ± 0.11 (effective charge = Z e jf) and 
5expUA = 3.02 ± 0.21 (effective charge = Z e jf — 1). 

We now argue why S exp Ua is indeed a relevant quantity, to a good approximation. At 
the 9 point we have B 2 = so that Jg = —8/3 from Eq. (J5J). From tables I and II, we 
see that generally \i ^> 1 and a( <C 1; hence, we have J\ ~ 4£//i and J 2 — 45 exp — £ for 
often fi5 > 1. This would lead to SexpUA — 4.4. On the other hand, at very high /, J\ 
and £ tend to zero and, because [/^ > 1, the scaled virial coefficient B 2 /B2 S reduces to 
— I J 2 exp [7,4 ~ —^5 expUA leading to 6expUA — 3 estimated from Fig. 3. Hence, the two 
estimates at the respective extremes are fairly consistent. To summarize, we may propose a 
crude approximation to the second virial coefficient which is a universal function of 6 exp Ua 

J^+I- V < 28 > 

The third term on the right is exact in the limit 5^-0, whereas the absolute error in the 
second term is smaller than 0.25 when J > 0.1 M. Using Eq. (|2*Hj) to fit the data leads to 
5 exp Ua = 4.2 when we use the effective charge Z e ff, whereas 5 exp Ua = 3.7 when we use 
the lower effective charge Z (see Fig. 5). 

In Fig. 3 we see that the curves at low values of 5 fit the data at high ionic strengths 
better. In the remainder of this article, we therefore employ the values 5 = 0.079 and 
Ua = 3.70, corresponding to the lowered effective charge Z and Ig = 0.21 M. In Fig. 6 we 
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show a comparison between the experimental data at a pH of about 7.5 and the theoretical 
curve computed numerically with the same parameters. 

3. AHS potential 

Values of eo, Us and r at several ionic strengths are given in tables I and II. Fig. 7 
displays the ionic-strength dependence of the adhesion parameter r. Near the 6 point, r 
decreases quickly with increasing I. At high ionic strength, r approaches the limiting value 
of (65{e UA - 1)) , which, upon the use of our choice 5 = 0.079 and Ua = 3.7, is equal to 
0.0535. We note that at pH 4.5 and at ionic strengths I = 0.05 M and I = 0.1 M, the 
computed values of eo, Us and r become nonsensical. In that case, the attractive potential 
is simply not strong enough to compensate the electrostatic repulsion completely so our 
analytical approach breaks down. This can also be seen in Fig. 2, where we have B 2 > B 2 S 
for these two values of the ionic strength. The same effect occurs at pH 7.5 when / = 0.05 
M. 



III. LIQUID STATE THEORY AT HIGHER DENSITIES 
A. Theory 

1. Density dependent attractive well in the Baxter limit 

In section II, we introduced the AHS potential as a convenient first approximation to the 
interaction between proteins. We determined the adhesion parameter r by matching values 
of the second virial coefficient which is methodologically correct only in the asymptotic limit 
of very low densities. In this section we propose a new procedure of choosing r, which is 
valid at higher densities but r now depends on the protein density. We extend a method 
originally proposed by Weeks, Chandler and Anderson j^J for repulsive interactions. They 
variationally determined an effective hard sphere diameter for a soft, repulsive potential of 
short-range, but we argue that their scheme is more generally applicable as long as the full 
interaction remains of short range, which is the case here. 

We start by introducing a functional expansion of the excess Helmholtz free energy AA 
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in terms of the Mayer function of the interaction U 

p~ x A(p, T; <p 8 ) = p-'Aip, T- cpahs) + | A J rfx B AHS {x)+ (29) 
i] 2 f 3 \ 2 a 3 f (3 

+ T V 47T j V / rfXlrfX2C?X3 ^^(s12)5aH5(^13)^ 3 H5( x 1 ) x 2,X 3 ) + . . . . 

Here, V is the volume of the system, .4. = — AA/V, f s (x) = e~ u ^ x \ ^>ahs( x ) — e~ UAHS ^ x \ 
rj = 47ra 3 p/3 is the volume fraction of particles, ./',V f ^ ( x i , x?, x a) i s a complicated function 
depending on two and three particle correlation functions (see |3(l ]) and x i2 = xi — x 2 etc. 
We define the quantity 

Bahs(x) = Vahs(x) [<p s {x) ~ <Pahs(x)} , (30) 

in terms of the so-called cavity function Hahs(x) = gAHs(x) / ^Pahs(x) = ^-^^y and radial 
distribution function gAHs{ x ) pertaining to an appropriate AHS potential which is the ref- 
erence state. Both these functions depend on p, T and the effective adhesive parameter 
r, the latter to be determined variationally. From now on, we omit the subscript AHS in 
B A hs( x ), 9ahs( x ) etc. for the sake of brevity. 

We next choose r by requiring that the first-order correction to the excess free energy 
vanishes 

/ dxB(x) = 0. (31) 

This is the analogue of Eq. ()12j) . Hence, in the spirit of the previous section, we split up 
this integral into two parts. The first indicates that the tail of the electrostatic interaction 
is compensated by part of the original square well attraction 



oo 

2 



dxx 2 B(x) = (32) 

'2+e 

(0 < e < 5) and yields e. The second determines the density dependent strength r of the 
AHS interaction 

2+e 

dx x 2 B(x) = 0. (33) 

This expresses the fact that the reference potential has to compensate for the remaining 
part of the original interaction. 
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2. Approximate radial distribution function for the Baxter potential 



In order to be able to determine r from Eqs. (}3*2"j) and (|33|). we need to know g(x), the 
radial distribution function of the reference interaction, the AHS potential. In the Percus- 
Yevick approximation developed by Baxter, g(x) has a singular contribution g w {x) which, 
after the limit uo — > 0, acts like a delta function and results from the stickiness of the 
interaction at the surfaces of two touching spheres. We split g(x) into g u (x) and a regular 
term g(x) 



with 



g{x) = g{x) +g w {x) 







X(2+u>) 
V2w 



x < 2 

+ 0(1) 2 < x < 2 + u) 
x > 2 + u 

analogously to Eq. (|24j). where the amplitude A is the smaller of the two solutions of 



(34) 



(35) 



T 



1+7//2 1 



n 



(36) 



;i-?7) 2 A 1-77 12' 

For x < 2, <?(;r) equals zero owing to the hard-core repulsion, whereas g(x) tends to unity 
for large x. For proteins, it turns out that <p s (x) — (Pahs(%) is often appreciably nonzero 
only near the surface of the sphere so we approximate g(x) in the interval 2 < x < 4 by the 
first two terms of its Taylor expansion 



9{x) 



x < 2 
G(l + H(x-2)) 2<x<4 

1 x > 4 



(37) 



The two constants have been derived by Bravo Yuste and Santos |31l 



G = Ar 



(38) 



and 



H 



V 



1 + 11*7, 1 
-A + 



577 9(1 + 77) 1 



;i-77) V 12 " 12 "'i-77 2(1 - 1]) 2 x 1 ' (o<)) 

Numerical work j32j bears out that Eq. ()37|) is quite reasonable since the range of both 
attractive and electrostatic forces is much smaller than the diameter of the protein. 
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3. Determination of the effective adhesion 



We next determine r from Eq. ()33|). first using Eq. (|32j) to obtain e. From Eqs. 
(I3(J|) and (J34|) . the function 5(a;) can be shown to have the following form (repressing terms 
that ultimately disappear in the limit uj — > 0) 

= 5(x)-^(x), (40) 

where the regular term is given by 

,0 < x < 2 

B(x) = { ~ ~ . (41) 

(g-tf(*) _ \)g{x) X >2 



poo r2+8 poo 

/ dxx 2 B(x) = I dxx 2 B(x) + / dxx 2 B(x) = 0. (42) 

</2+e </2+<e J2+S 



Eq. (J3~2~j) may be conveniently expressed as 

t>2+& poo 

dx x 2 B(x) + I 

>2+e J 2+5 

Using dx f(x) ~ ~(5 — e) [f{2 + 5) + f(2 + e)] and neglecting terms of order 5 2 and e 2 , 
we write the first integral as 



PZ + O 

/ dxx 2 B{x) ~G{8- e)Kx{S,e), (43) 

J2+E 



with 



K x (8, e) = 2 (e^e-^ 6 ^ - 1 ) (1 + (1 + H) 5) + 2 (e^e - ^ 6 - l) (1 + (1 + H) e) . 

(44) 

Again, we stress that, although 5 <C 1 and e « 1, |ii and fie may be of order unity. 
Furthermore, we note that if we take the limit rj J, 0, then A — > r _1 and Cr — > 1, so we 
recover Eq. (|16|) if we neglect terms of order 5 and e. We tackle the second integral by 
adopting the approximation: 1 — exp(— U(x)) = 1 — exp(2^x~ 1 e^^ x ^) ~ 2^x^ 1 e^^ x ~ 2 ^ — 
2£ 2 x~ 2 e~ 2 ^ x ~ 2 ' > + 2£ 3 x~ 2 e~ 3fJ "( x ~ 2 ^ /3 (note that in this Taylor expansion of the exponential 
we have replaced one factor x^ 1 by 2~ l in the last term). We then write 

POO 

- dxx 2 B(x)~G((l + 5H)P 1 + HP 2 ) (45) 

J2+5 



with 



f°° 8 16 / 8 \ 

Pi=/ dxx 2 (l-e- u{x) ) ~ — (1 + //5)M + — M 1 - M + -M 2 (46) 
J2+5 ^ ^ V 9 / 
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and 

roc o 1 fi / 1 R \ 

P 2 = / dx x 2 (x — 2 — <5)(1 — e~ u ^) ~ —(2 + fi5)M H — - ( M — -M 2 + — M 3 ) . (47) 

m V 2 27 y 

Here, M = fe _/ "74. Using the approximations 1-M + 8M/9 ~ (1 + M) -1 and M-M 2 /2 + 
8M 3 /27 ~ log(l + M) , we arrive at 

p ^> + ^ m + 7TTm < 48 > 

and 

P 2 ~ A(2 + ^)M + ^log(l + M). (49) 

Hence, the variable e, which depends on the density by virtue of the density dependence of 
G and H, is determined iteratively from 

(l + ^)P 1 + gP 2 
(5 - e new = . (50) 

One starts with e ^ = 5 and iterates until a stationary e new is reached. 

The next step is to calculate r from Eq. (JEBJ), which, with the help of Eq. (|4*0*j) . is 
equivalent to the expression 

dxx 2 B{x) = —. (51) 



We have taken the limit u — > 0. Again using the approximation J^ +e dx /( 



|e[/(2 + e)+/(2)], we write 

/■2+e _ _ 

y d« 2 5(i)~2Ge [^e"^ - lj (1 + (1 + #) e) + (e^e - * - l)J . (52) 
Together with the expressions Eq. (JoTj) and G = At (Eq. (|35jl ). this leads to 

(>A e -Trb e ~ Me - l) (1 + (1 + H) e) + (e UA e-* - l)] . (53) 



1 

— ~ 3e 

T 



Accordingly, r may be determined iteratively if we recall that both H and e also depend 
on r. A way of quickly determining r and e is choosing a starting value for both (e = 5 
and r = 0.2 say), and then alternately using Eqs. (f5U)) and (|53~|) until the iterates become 
stationary. 
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B. Application to lysozyme 



We have already determined the interaction in section IIB2b. We next compute r itera- 
tively and it now depends on both the density of protein and the ionic strength. (See table 
III). 

Thermodynamic properties like the osmotic compressibility kt are also simply obtained 
from r. In the Percus-Yevick approximation, kt is given by [19] 

(ok Tk r 1 - 1 9U - (l + 2q-Aq(l-q)) 2 

where A is the smaller of the two solutions of Eq. (jHfijl . Fig. 8 compares the predicted 
density dependence of the (scaled) inverse osmotic compressibility at various ionic strengths 
with experimental data from H and Q- 

IV. DISCUSSION 

One difficulty in comparing our computations with experiment has been the substantial 
margin of error in the osmotic measurements. In the case of other biomacromolecules like 
rodlike DNA, it has been possible to obtain the second virial B 2 at better than 10% accuracy 
H H Q. One possibility for the occurrence of discrepancies in B 2 is the variety of 
lysozyme types. Poznanski et al 36J have established that popular commercial lysozyme 
preparations like Seikagaku and Sigma exhibit significant differences under dynamic light 
scattering. Nevertheless, the variation in B 2 at, say, about 0.5 M NaCl (see Fig. 3), is so 
large that it needs to be explained. 

The relatively large variation in the experimental measurements of B 2 makes it difficult 
to falsify stringently other models of attractive forces like that of van der Waals type, for in- 
stance. It proves feasible to get satisfactory agreement with the experimental data displayed 
in Fig. 3 if we let the dispersion interaction be given by the nonretarded Hamaker poten- 
tial [18] for spheres of dimensions appropriate for lysozyme, with an adjustable Hamaker 
constant of order k^T though with a very short cut-off at around 0.1 — 0.2 nm. However, 
the necessity of such a cut-off, which is already beyond the limit of validity of continuum 
approximations, may be viewed as positing the equivalent of a short-range interaction like 
that of Eq. (jHJ), in large part. The long-range dispersion interaction beyond some distance 
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much smaller than the radius a, plays only a minor role. 

Ste n has c^ed^eB^^ because dive, gences in the free energy appear at 
the level of the 12th virial. Therefore, the most straightforward way to interpret our liquid 
state theory is to stress that our zero-order theory describes the reference state only up to 
and including the 11th virial within the Percus-Yevick approximation. The analysis of phase 
transitions must be viewed with caution (for a comparison of recent simulations — taking 
the limit of zero polydispersity after the limit of vanishing well depth — with Percus-Yevick 

n 

theory, see [38]). A second problem is here that, at large ionic strengths, a considerable 
electrostatic repulsion is balanced against a significant attraction (see Fig. 1) and it is 
difficult to see how good such a compensatory scheme should work at high concentrations 
near dense packing. 

In summary, we have presented a fairly good theory of the ionic-strength dependence of 
the osmotic properties of lysozyme in terms of a sticky interaction which is independent of 
charge or salt concentration. This conclusion, by itself, is not new for it has been reached 
earlier by formulating numerical work incorporating short-range forces and screened elec- 
rostatics and comparing it with X-ray scattering [3, 13 an d liquid-liquid phase separation 
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merit of the current analysis is its transparency because it is analytical 
and it is based on a nonperturbative variational principle for general short-range potentials 
so it may be readily generalized. 



V. APPENDIX 

A. Effective charge 

For the repulsive tail of the two particle interaction, we use the Debye-Hiickel poten- 
tial, which is the far-field solution of the Poisson-Boltzmann equation. In our case, the 
(dimensionless) potential at the surface is often merely of order unity, so the Debye-Hiickel 
potential slightly overestimates the solution to the Poisson-Boltzmann equation. To remedy 
this, we use a renormalized charge within the Debye-Hiickel potential, chosen in such a way 
that, at large distances, the Debye-Hiickel potential coincides with the tail of the solution 
of the Poisson-Boltzmann equation determined by the real charge 0|. This will result in 
an underestimation of the potential at small separations, but the form of the Debye-Hiickel 
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potential we use here (Eq. is in fact only accurate at large separations and overestimates 
the interaction at small separations appreciably i.e. when overlap of the two double layers 
occurs (by about 20%, see [l8|). The two effects thus partly cancel, although the latter 
effect is larger than the former. 

The Poisson-Boltzmann equation for the dimensionless potential if) (r) = q<f)(r) jk^T of a 
single sphere of radius a and total charge qZ, assumed positive for convenience, immersed 
in a solvent with Bjerrum length Q, at a concentration of ions leading to a Debye length k, 
is written as 



— — r 2 — if) (r) = k 2 sinh if) (r) , 
r 2 dr dr 



with boundary conditions 



^-if) (r) 
dr 



— ^; lim ip(r) = 0. 



(55) 



(56) 



Linearizing Eq. (}55|) (^) < 1), we find the Debye-Htickel solution 

ZO e - K ( r - Q ) 

= • (57) 

1 + n r 

We next derive the first-order correction to this solution. Putting if) (r) = ipo (r) + if>x (r), 
with \ipi (r)| <C IV'o ( r )|> results in the following linear differential equation for if)\ 



\ d n d 1 q / S / •! 

^dT ^ (r) = 6" ^ (r) • 



(58) 



Keeping in mind that if)\ (r) = o(^o ( r )) ; as r ^ oo, we integrate the differential equation 
once to obtain 



d 
dr 



ipx (r) 



k 2 ( ZQe^y Ei(3nr 



6 + 



and a second time to derive 



, . . k 3 ( ZQe» 
Wi \r) = — — 



*(— -f3 + -)£ l (3/ 1 ;; 
6\l + yU/\/tr \ nr 



(59) 



(60) 



where -Ei(x) is the exponential integral defined by E\(x) = dtt 1 e t . Using the first of 
the two boundary conditions, we then determine the renormalized charge Z e ff 



7 ° 2 d i n 

= z-JL(9. 



r=a 
2 



: Q ~dr 

Z 
1 + n 



"00 (r) 



a 2 d , , . 



(61) 
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where 

F(jm) = 3/ie^(3A*) ~ 1 - ±- + ^ - . . . (62) 

Recapitulating, we have calculated, to leading order, the charge Z e ff which has to be 
inserted into the Debye-Hiickel potential (Eq. (0])) so that this has the correct asymptotic 
behavior at large r, coinciding with the tail of the Poisson-Boltzmann solution. 



B. Dependence of B 2 on ionic strength 

Here, we prove some simple inequalities describing the behavior of the second virial 
coefficient as a function of the ionic strength for an interaction consisting of a Debye-Hiickel 
repulsion Udh{x) and a general attractive potential Ua(x), the latter not depending on the 
ionic strength. If we let U{x) = Udh(x) + Ua(x), then B 2 is given by Eq. © with 



J = J dxx 2 (l- e- u &) . (63) 

Then, we have 

d 4- = r dx , 2 ^M e -^) = r dx x 2 ( ^ - (x - 2)) U DH (x) e-^). (64) 
dfi J 2 d/i J 2 \ dfi J 

In Fig. 9 we see that in the regime of interest < 0, so we conclude that 

dB* = 3 B HsdJ < Q (65) 
dfi 8 d\i 

In the same way it is clear from the second derivative 
d 2 J f°° 2 {d 2 \n£ fdhxi x2 



dfi 2 



r dx x2 ( + ■ {x ~ 2) ) (i _ udh {x)) ) udh {x) e ~ u(x) (66) 



and the fact that > in the regime of interest, that 



d/i 2 8 d/i 2 
if Udh (2) < 1, i.e. if £ < 1 (a sufficient condition). 



- B* s —>0, (67) 



C. Corrections to the free energy 

In section III, we viewed a suspension of proteins as a system of spheres with an AHS 
interaction and we chose the parameter r of the AHS potential such that the first order 
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correction in the functional expansion of the free energy (Eq. (|29j)) vanishes (see Eq. ()3ip). 
In an attempt to justify this approximation and explore its regime of applicability, we 
estimate the size of the second order correction to the free energy (from Eq. (|29|) ) which is 
either positive or negative definite 



A = i aV" 1 J dxndx 2 dx 3 B(x 12 )B 



9 

(xi 3 )h(x 23 ) = jrfY. 



(68) 



It is convenient to rewrite the integral in such a way that the angular integration can be 
performed explicitly (see below). 



/*0O /*0O /*7T 

Y= dtt 2 B(t) / dss 2 B(s) / d$ sin ■& h(Vs 2 + t 2 - 2st cos tf) 
Jo Jo Jo 



(69) 



s+t 



= 2 1 dttB{t) I ds sB(s) / duuh(u). 

Jo Jt Js-t 

Here we have used the Kirkwood superposition approximation J^ 3 ] f ( X i, x 2 , x 3 ) = h(x 23 ) [30], 
where h(x) = g(x) — 1 is the pair correlation function. We have employed the substitution 
u 2 = s 2 + 1 2 — 2st cos with $ the angle between x 12 and x 13 . Using the expression for g(x) 
(Eq. (|34p) and defining /i(x) = #(a;) — 1, we split Y into three parts 

Y = Y + Y 1 + Y 2 (70) 

where we have introduced the limit u — > and where 

-t+2 



Y = — dt tB{t) / rfs s£(s) ~ 

3 J 2 A 

o \ roc roc \ |~ /*oo 

~ — / rft t£(t) / ds sB(s) = - dt tB(t) 
3 J 2 Jt 3 \_J 2 



(71) 



and 



/>oo rt+2 rs+t 

Yi = 2 / dttB(i) / dssB(s) / duuh(u). 

J2 Jt Js-t 

roc rt+A fs+t 

Y 2 = 2 dt tB(t) / s5(s) / dw u/i(u) < K a 

</2 A+2 A-t 



(72) 



(73) 



To simplify Eq. ((721), we substitute Eq. fl37J) and note that s + t>4and0<s-t<2. We 
then derive 

du uh{u) = ~ (9G + 10G# - 12) + hs - t) 2 . (74) 
s-t 3 2 
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Next, using Eq. (|3ip. we integrate the nonconstant term leading to a product of two integrals 



/oo rt+2 poo poo 

dttB(t) ds sB(s)(s-t) 2 ~ y dttB(t) y ds sE 



(s)(s-tf 



(75) 



dt t5(t) 



ds s 3 B(s) 



Hence, Y\ is written in terms of one- dimensional integrals 

1 2 



Yi~- (9G + 10GH - 12) 



and this is also the case for Y 



Y~ - ( 9G+ lOGTd- 12 + - 
3 V 2 



dt tB(t) 



dt tB(t) 



dt tB(t) 



+ 



dt tB(t) 



ds s 3 B(s) 



/"OO 

J ds s 3 B(s) 



(76) 



(77) 



Our goal is to obtain explicit approximations for these integrals by expediently using Eqs. 
(|32j) and (jHEJ). First, we consider integrals on the interval [2, 2 + e] which are dominated by 
the singular part of B(x). We substitute Eq. (|4T)|) into (J3*H|) and let lj — ► 



2+e ~ 9 X 
dt t 2 B(t) = — . 
2 3 



(78) 



We use this relation to rewrite part of one of the integrals in Eq. (|77|) in two ways, noting 
that e < 1. 



2 + 6 



2+e 



2+e 



2+e 



dttB(t) = --+l dttB{t) = --j dt {t - 2)tB(t) ~ -- / dttB(t). (79) 



We thus conclude that 



2+e 



dt £B(t) 



4/ 3 



(80) 



so the first equality in Eq. (J79*j) allows us to attain the explicit expression 



/•2+e 

y dt t.B(t) ~ - 



Ae 
12' 



(81) 



Similarly, we use Eqs. (jjUj) and (|75j) to evaluate part of the other integral in Eq. (|77|) . 

2+e A \ /-2+e /-2+e \ 



dt t d B(t) = + 



dt t 3 B(t) 



2 3 J 2 </2 3 



dt (t-2)t 2 B(t) ~ — . 



(82) 



We note that both integrals in Eqs. (|8*T|) and (jBZj) are 0(e) because the integral in Eq. (|7H|) 
is independent of e owing to the singular part of B(x). If -B(x) had been completely regular, 
the integrals in Eqs. (}8*T]) and (|52*|) would have been 0(e 2 ). 
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We next consider the remaining two integrals on the interval [2 + e, oo). We start by 
splitting Eq. (}3*2*)) into two parts since 2 + 5 demarcates two different regimes 

p2+8 poo 

/ dt t 2 B(t) + dt t 2 B(t) = 0. (83) 

J2+e J2+8 

Using this equation and the approximation B{t) ~ — 2£e _M (*~ 2 < ) /t, we may simplify the two 
integrals, ultimately omitting 0(5) terms 

-2+8 



poo pi+o poo 

j dt tB(t) = dt tB(t) + dt tB(t) 

J2+e J2+e J2+8 



1 A 5 



2 V 2 



2+8 



1 - - / dt t 2 B(t) + / dt tB{t) 



2+e J 2+8 

oo i roo 



- \ dt t 2 B(t) - - / dttit- 2)B(t) ~ -^e"^ (84) 

4 J 2+8 2 J 2+s jj, 1 



poo p2-\-S poo 

/ dt t 3 B(t) = dt t 3 B(t) + dt i 3 B(t) 

J2+e J2+e J 2+8 

-2+8 



pl+d poo 

2 + 5) dtt 2 B(t) + / dtt 3 B(t) 

J2+e J2+8 



POO POO £ 

-5 dtt 2 B(t)+ dtt 2 {t-2)B{t)~-A\e~^. (85) 

J2+8 J2+8 f< 



We remark that both expressions in Eqs. and (J85j) are 0(fi 2 ) because B(t) is regular 
for t > 2 + e. We then combine Eqs. (JHU) and JED, and Eqs. JH2J and JH^J 



(it £B(t) ~ + -^e-" 5 — — — / ds s 3 B(s). (86) 



2 12 /i 2 4 A 

Finally, using Eqs. (|B5j) . (J77|) and (J85J) . we arrive at an approximation for the correction to 
the free energy 

(87) 



A = \ 2 Y ~l v 2 (g + h-6+- 



4 ' 2 ' V 6 



/i 2 12 



Despite the variety of approximations used, this expression still retains its "definite" charac- 
ter (it turns out to be negative in the numerical calculations below). However, the numerical 
coefficients within the last quadratic factor are not exact. Furthermore, the status of the 
present theory differs from that of the Weeks-Chandler- Anderson theory j^J. In the latter, 
A is of fourth order in the perturbation whereas it is basically quadratic here for the reason 
stated below Eq. (J32J) . 
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To estimate the importance of this correction, we first calculate the osmotic pressure 
resulting from the neglect of second and higher order terms in the functional expansion Eq. 
(J2HJ). This amounts to determining r from Eqs. (J35J) . (|5L)|) and (J55|) and then computing the 
osmotic pressure from Ref. Q| 

n l+r] + r] 2 -\r](l-r])(l + ^r]) + X 3 r] 2 (l-r]) 3 /36 

. (88) 



pk B T (1 - r]) 3 

Then, we evaluate the correction to the osmotic pressure due to the second order term in 
Eq. ()29|). The osmotic pressure is related to the free energy by 

11 -, d -K^. (89) 



pksT drj 

Because Y depends only weakly on rj, we approximate the correction to the osmotic pressure 
by 

— ti— — ~ -2A. 90 
drj 

We have compiled the pressure and its correction in table IV for the same sets of parameters 
as in table I (omitting the trivial case where r] = 0). 
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Tables 
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l.Z I 


1 m 
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n HQ A 


U. f oz 


n /i no 
U.4U9 


U.Z9D 


U.ZZ9 


n 


1.25 


1.76 


2.16 


2.50 


2.79 


3.06 


3.74 


5.58 


6.83 


7.89 


eo 






0.0208 


0.0466 


0.0585 


0.0644 


0.0720 


0.0773 


0.0782 


0.0785 


Us 






2.26 


2.52 


2.70 


2.82 


3.05 


3.37 


3.47 


3.53 


T 






0.933 


0.314 


0.205 


0.164 


0.115 


0.0767 


0.0684 


0.0642 



Table I: Values of the actual charge Z (from [29]), the renormalized or effective charge 
Z e ff (from Eq. ()61)l ). the lowered effective charge Z = Z e ff — 1, and dimensionless interaction 
parameters £ and /z, and eo, Us and r as a function of the ionic strength /. The pH equals 
4.5 and £ has been calculated using the lowered effective charge Z. Values of Us and r have 
been computed using Eqs. (J2~3j) and (}2*T|) . respectively, and eo has been calculated using the 
procedure described immediately after Eq. (JTTj). 
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0.0889 


u 


1.25 
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0.0764 


0.0784 


0.0787 


0.0788 


U s 




2.83 


3.03 


3.14 


3.23 


3.28 


3.39 


3.56 


3.61 


3.63 


T 




0.212 


0.132 


0.108 


0.0943 


0.0877 


0.0758 


0.0623 


0.0590 


0.0574 



Table II: Same as table I, but now with a pH equal to 7.5. 
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0.15 M 


0.2 M 


0.25 M 


0.3 M 


0.45 M 


1 M 


1.5 M 


2 M 


o 


T 


0.829 


0.295 


0.194 


0.156 


0.110 


0.0735 


0.0656 


0.0616 




£ 


0230 

V J a V 7 aU '.J V 7 


0.0483 


0596 

V 7 • V 7 7 '7 V 7 


0653 

V 7 • V 7 V7 .7 '.7 


0.0725 


0.0775 


0.0782 


0.0786 


0.05 




0.712 


0.289 


0.193 


0.155 


0.110 










£ 


0.0266 


0.0492 


0.0600 


0.0655 


0.0725 








0.1 


1 


0.620 


0.283 


0.192 


0.155 


0.110 










£ 


0.0303 


0.0502 


0.0603 


0.0656 


0.0725 








0.15 


T 


0.544 


0.276 


0.191 


0.155 


0.110 










£ 


0.0342 


0.0514 


0.0607 


0.0657 


0.0724 








0.2 


J 


0.482 


0.268 


0.190 


0.155 


0.110 










£ 


0383 


0528 

V 7 • V 7 '.7 v7 


0.0611 


0658 

V 7 • V 7 V7 ' 7 v.7 


0.0723 








0.3 


T 


0.380 


0.251 


0.186 


0.154 


0.110 










e 


0.0477 


0.0563 


0.0624 


0.0663 


0.0722 








0.4 


T 


0.300 


0.228 


0.179 


0.152 


0.110 










e 


0.0600 


0.0619 


0.0651 


0.0677 


0.0724 









Table III: The scaled depth e of the effective attractive well and the strength of the 
effective adhesive interaction r at pH 4.5 as a function of the ionic strength / and volume 
fraction 77. The values of e and r have been evaluated from Eqs. (|5Uj) and ([5l5|) 
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0.15 M 


0.2 M 


0.25 M 


0.3 M 


0.45 M 


0.05 


n 

pk B T 


1.143 


1.033 


0.949 


0.889 


0.763 




-2A 


0.004 


0.001 


0.0004 


0.0001 


000008 

V J a V 7 V J V 7 V 7 V 7 V7 


0.1 


n 

pk B T 


1.290 


1.074 


0.915 


0.805 


0.575 




-2A 


0.019 


0.006 


0.002 


0.0005 


0.00002 


0.15 


n 

pk B T 


1.437 


1.123 


0.898 


0.749 


0.448 




-2A 


0.044 


0.014 


0.004 


0.001 


0.00002 


0.2 


n 

pk B T 


1.583 


1.183 


0.904 


0.721 


0.375 




-2A 


0.085 


0.026 


0.008 


0.003 


000007 

V 7 • V 7 V 7 V 7 V 7 V 7 1 


0.3 


n 

pk B T 


1.866 


1.361 


0.988 


0.753 


0.340 




-2A 


0.228 


0.068 


0.022 


0.008 


0.00002 


0.4 


n 

pk B T 


2.17 


1.659 


1.231 


0.960 


0.470 




-2A 


0.488 


0.143 


0.046 


0.016 


0.0001 



Table IV: The osmotic pressure from Eq. (J%B*j) and its correction from Eq. (190(1 as a 
function of the ionic strength I and the packing fraction 77. 
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Figure Captions 



Fig. 1: The integrand of Eq. (j2J) versus the distance r. As shown by the shaded regions, 
the repulsive tail is compensated by part of the attractive interaction. 

Fig. 2: Experimental data of the second virial coefficient B 2 of lysozyme as a function of 
the ionic strength J at a pH of about 4.5. The second virial coefficient is scaled by the hard 



sphere value B 2 S ■ Black squares: Bonnete et al. |27[ 



et al. 
et al. 



23 1, pH 4.5, 20°C; grey squares: Muschol et al. 



pH 4.5, 20°C; grey triangles: Curtis 



24j^_pH 4.7, 20°C; black stars: Curtis 
2J, pH 4.5, 25°C; black triangles: 
Velev et al. |21[, pH 4.5, 25°C; white squares: Rosenbaum et al. 2(|, pH 4.6, 25°C; white 



22], pH 4.5, 25°C; black diamonds: Bonnete et al. 



diamonds: Rosenbaum et al 
25°C; white stars: Piazza et al. 



14|L_P H 4.6, 25°C; grey stars: Bloustine et al. Q>£H 4.6, 
2J,_pH 4.7, 25°C; white triangles: Behlke et al. [28], pH 



4.5; grey diamonds: Bloustine et al. [26j, pH 4.7. In all cases, the electrolyte is NaCl, often 
with a small amount of Na acetate added. 

Fig. 3: A fit of Eq. El to the experimental data of Fig. 2 (except for those of Refs. 
and 28]). On the right-hand side of the figure, the upper solid line corresponds to Ig = 0.19, 
5 = 0.564 and U A = 1.48, the upper dotted line to Ig = 0.20, 5 = 0.468 and U A = 1.70 and 
the middle solid line to Ig = 0.21, 5 = 0.379 and Ua = 1.95, all at an effective charge Z e ff. 
The middle dotted line corresponds to Ig = 0.19, 5 = 0.25 and Ua = 2.4, the lower solid one 
to Ig = 0.20, 5 = 0.167 and U A = 2.87 and the lower dotted one to Ig = 0.21, 5 = 0.079 and 
Ua = 3.70, all at a lowered effective charge Z. 

Fig. 4: Experimental data of the second virial coefficient B 2 of lysozyme as a function 
of the ionic strength / at a pH of about 7.5. The second virial coefficient is scaled by the 
hard sphere value B^ 3 . Black stars: Rosenbaum et al. [2^, pH 7.4, 25°C; black triangles: 
Velev et al. |2l|, pH 7.5, 25°C; black squares: Rosenbaum et al. [14J, pH 7.8, 25°C. 

Fig. 5: Fits of Eq. (|2*H|) to experimental data of Fig. 3. Full line {Z e ff and SexpUA = 
4.2); Dotted line (Z and 5exp U A = 3.7). 

Fig. 6: Comparison between the experimental data at pH 7.5 and full theory Eq. 
Parameters as in the lower dotted curve in Fig. 3 (5 = 0.079 and Ua = 3.70). 

Fig. 7: Ionic-strength dependence of AHS parameter r at pH 4.5 and pH 7.5. The 
dotted line denotes the limiting value of r as / — » 00. 

Fig. 8: Inverse osmotic compressibility as a function of the volume fraction rj at various 
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ionic strengths. Experimental data: black squares: / = 0.18 M; black triangles: I = 0.23 
M; black stars: / = 0.28 M; black diamonds: / = 0.33 M; open squares: / = 0.48 M. All 
data from Rosenbaum et al. l^J, except for those at I = 0.23 M (black triangles) (Piazza 
et al. ^J). Curves computed from Eq. (J54j) with 5 = 0.079, U a = 3.70 and the lowered 
effective charge Z; r has been determined from Eq. (|53jl . From top to bottom: / = 0.18 M, 
/ = 0.23 M, I = 0.28 M, I = 0.33 M and / = 0.48 M. 

Fig. 9: Dependence of ln£ on fi at pH 4.5 and pH 7.5. In both cases dhx^/d^i < and 
d 2 ln£/d/i 2 > if 1 < fi < 8. 
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